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A TIGHT Liso - Lp.z - T 0 CORRELATION OF GAMMA-RAY BURSTS 
En-Wei Liang 1,2,3 , Ting-Ting Lin 1,2 , Jing Li) 1,2 , Rui-Jing Lu 1,2 , Jin Zhang 3,2 , Bing Zhang 4,1,2 

ABSTRACT 

We select a sample of 34 gamma-ray bursts (GRBs) whose T 0 values are derived with the onset 
peaks observed in the afterglow lightcurves (except for GRB 060218 whose To is estimated with its 
radio data), and investigate the correlations among To, the isotropic peak luminosity (Lj so ), and the 
peak energy (L PjZ ) of the vf v spectrum in the cosmological rest frame. An analysis of pair correlations 
among these observables well confirms the results reported by the previous papers. More interestingly, 
a tight correlation among Li so , E P:Z , and To is found from a multiple regression analysis, which takes 
the form of L iso oc £i-34±o.i4 r i.32±o.i9 Qr ^ L 0.55±0.06p-0.50±0.u Nine other GRBs whose ro 

are derived via the pair production opacity constraint also follow such a correlation. Excluding GRB 
060218, the L; so — E PiZ — To correlation is valid, and it even holds in the jet co-moving frame. However, 
GRB 060218 deviates the L iso — E p relation of typical GRBs in the jet co-moving frame with 3cr. We 
argue that the Li so — E P}Z — To correlation may be more physical than the L lso — L PjZ correlation, since 
physically the relationship between the observed L ISO and E PyZ not only depends on radiation physics, 
but also depends on the bulk motion of the jet. We explore the possible origins of this correlation and 
discuss its physical implications for understanding GRB jet composition and radiation mechanism. 
Subject headings: gamma-ray burst: general-methods: statistical 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous elec¬ 
tromagnetic explosions in the Universe. The observed 
GRB spectra are typically fit by an empirical smoothly- 
jointed broken power-law function, the so-called Band 
function, characterized by a peak energy (E p ) in the vf v 
spectrum, which ranges from keVs to MeVs (Band et al. 
1993; Zhang et al. 2011). The physical origin of the E p 
and the Band function spectrum is still subject to de¬ 
bate (Kumar & Zhang 2015 for a recent review). One 
possibility is synchrotron radiation in internal shocks 
(e.g., Meszaros et al. 1994; Daigne & Mochkovitch 1998; 
Daigne et al. 2011) or in internal magnetic dissipation 
regions (e.g., Zhang & Yan 2011). A recent development 
of this model suggests that if the emission region is far 
from the central engine, the predicted spectra mimic the 
observed Band spectra (Uhm & Zhang 2014), and can 
well fit the data with comparable confidence level as the 
Band function (Zhang et al. 2015). Another scenario 
attributes the E p and Band function to emission from 
a dissipative photosphere of the fireball (e.g., Rees & 
Meszaros 2005; Beloborodov 2010; Lazzati & Begelman 
2010; Lundman et al. 2013). These models invoke dif¬ 
ferent compositions of GRB jets (e.g., matter dominated 
fireballs vs. Poynting-flux dominated jets). Differenti¬ 
ating them from the data has profound implications in 
understanding the physical mechanisms of GRBs. 

Several important empirical correlations have been dis¬ 
covered with the GRB data. These correlations not only 
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give important clues to understand GRB physics (even 
though some correlations still lack straightforward theo¬ 
retical interpretations), but also some of them have been 
used to constrain cosmological parameters (e.g., Wang 
et al. 2015 for a recent review). For example, Am- 
ati et al. (2002) found a tight correlation between the 
isotropic gamma-ray energy Ei so and the cosmological 
rest-frame peak energy E PtZ = E p ( 1 + z ). A similar cor¬ 
relation was found between the isotropic luminosity and 
E PtZ , both among different bursts (Yonetoku et al. 2004) 
and within a same burst (Liang et al. 2004). Ghirlanda 
et al. (2004) found a tighter relation by replacing E[ so 
with the geometrically-corrected energy of the GRB jets. 
The initial Lorentz factor (To) is a crucial parameter to 
understand GRB physics. Liang et al. (2010) discovered 
a correlation between (T 0 ) and U; so (see also Ghirlanda 
et al. 2012). Lii et al. (2012) showed that the isotropic 
luminosity (Li so ) is also correlated with To- 

Theoretically, the predicted E p , z not only depends on 
luminosity Li so , but also depends on the bulk Lorentz 
factor T 0 of the outflow (see e.g., Table 1 of Zhang & 
Meszaros 2002 for a summary). Therefore, it is inter¬ 
esting to search for possible multi-variable correlation 
among L PiZ , L lso (or Ei so ) and To- Some multi-variable 
correlations have been found for GRBs (e.g., Liang & 
Zhang 2005; Rossi et al. 2008), which are useful to un¬ 
derstand GRB physics. 

This paper presents a multiple-variable regression anal¬ 
ysis among E\ so (or L- lso ), E P}Z , and To for long duration 
GRBs. We compile a sample of long GRBs whose E- lso (or 
Lj so ), E p ^ z , and To can be derived from the observation 
data (Section 2). Several correlations among these quan¬ 
tities are presented in Section 3. Physical implications 
are discussed in Section 4, and conclusions are drawn in 
Section 5. Notation Q n = Q/ 10” is adopted in cgs units. 

2. DATA 
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We compile a sample of GRBs whose E- lso , Lj so , E PiZ , 
and To are available in the literature or can be calcu¬ 
lated with observational data. E p>z can be obtained by 
the measured E p and redshift 2. Both E \ so and L- lso are 
corrected to an energy band of 1 — 10 4 keV in the burst 
rest frame. 

The Lorentz factor To is a key parameter in this anal¬ 
ysis. Three methods have been proposed to estimate 
To of a GRB fireball. The first one is to use a smooth 
onset bump observed in optical afterglow lightcurves. 
By interpreting this bump as a result of deceleration 
of the fireball by an ambient medium in the thin shell 
regime (which is usually satisfied for a constant density 
medium), one can estimate To with the peak time of the 
bump (Meszaros & Rees 1993; Sari & Piran 1999; Zhang 
et al. 2003; Molinari et al. 2007). The second method 
is based on the “compactness” argument by interpreting 
the high energy cutoff of the prompt gamma-ray spec¬ 
trum as pair production signature (e.g., Baring & Hard¬ 
ing 1997; Lithwick & Sari 2001; Gupta & Zhang 2008). 
The third approach is to use a blackbody component de¬ 
tected in the spectra of some GRBs (e.g., Pe’er et al. 
2007; Peng et al. 2014; Zou et al. 2015). This method is 
based on the assumption of a matter-dominated fireball, 
and therefore is not reliable if the central engine carries a 
significant fraction of Poynting flux (Gao & Zhang 2015). 

The first method gives the most robust estimates to 
Po, since the deceleration time weaklier depends on other 
model parameters than To- Such afterglow onset feature 
is observed in about 1/3 of GRBs with early optical after¬ 
glow observations (Li et al. 2012). Although the “onset 
afterglow” feature in some X-ray and GeV lightcurves is 
also observed (Xue et al. 2009; Ghisellini et al. 2010), 
early X-ray emission and GeV emission may have con¬ 
tamination from internal emission components. Thus, we 
only include GRBs with an onset bump detected in the 
optical band (most bursts from Liang et al. 2010, 2013 
and references therein). We calculate their To values with 
the afterglow onset peaks assuming that the fireballs of 
these GRBs are in an inter stellar medium with density 
profile n = 1 cm -3 and their radiation efficiencies are 
r] = 0.2 (e.g., Liang et al. 2013). The only exception 
is the nearby low luminosity GRB 060218 whose To is 
robustly estimated from the radio data (Soderberg et al. 
2006). The inclusion of this GRB stretches the dynami¬ 
cal range of the correlations significantly. L- lso values are 
measured in the Is peak time of their lightcurves, and 
E iso values are calculated with their gamma-ray fluences. 
Both I/i so and Ei so are corrected to 1 — 10 4 keV band in 
the burst frame. Their E p values are derived from the 
fits to their time-integrated spectra with the Band func- 
tior0 Altogether our sample includes 34 GRBs, which 
are reported in Table 1. 

3. CORRELATIONS AMONG Ei SO , Gso, Ep.z, AND T 0 

With the available data, we conduct both single and 
multiple variable regression analysis to look for corre¬ 
lations among the parameters E- lso , L- lso , E PjZ , and To- 
Notice that the results depend on the specification of de¬ 
pendent and independent variables (Isobe et al. 1990). 

5 The E p of time-integrated spectrum is usually consistent with 
that derived from the spectrum around the peak time (Lyu et al. 
2015). 


Therefore, one may find discrepancy of the dependences 
among variables by specifying different dependent vari¬ 
ables for a given data set, especially when the data have 
large error bars or large scatter. To avoid specifying in¬ 
dependent and dependent variables in the best linear fits, 
in principle the algorithm of the bisector of two ordinary 
least-squares may be adopted. However, in connecting 
physical models, often some parameters (e.g., y ) are be¬ 
lieved to depend on other parameters (e.g., x), see more 
discussion in Section 4. We therefore use the Spearman 
correlation analysis to search for pair correlations among 
these parameters, and adopt the stepwise regression anal¬ 
ysis method to perform a multiple regression analysis for 
multiple parameters by specifying a given y. We measure 
the dispersion (a) of a regression model with standard de¬ 
viation of y r from y 1 where r marks the y value derived 
from the regression model. 

Figure 1 shows the pair correlations among E- lso , E PiZ , 
and To, or among L[ so , E PjZ , and To- They are derived 
in the logarithmic space, and the results are reported in 
Table 2. The regression lines together with their 2<r dis¬ 
persion regions are also shown in Figure 1. Correlations 
among £) so , Li so , E PjZ , and Fo are found, with a Spear¬ 
man correlation coefficient p > 0.6 and chance probabil¬ 
ity p < 1CT 4 . Tight -Eiso - Ep, Eiso - E p , E iso - F 0 , 
E Pi is 0 — To correlations found in previous papers (e.g., 
Amati et al. 2002; Liang et al. 2004; Yonetoku et 
al. 2004; Liang et al. 2010; Lu et al. 2012; Lii et al. 
2012)0 are well confirmed. Combining the Ei so — E p z 
(or Eiso - E PiZ ) and the E iso - r 0 (or E iso - T 0 ) corre¬ 
lations, one may suspect a correlation between E PjZ and 
Fo- Such a correlation is indeed found in our sample, 
which is E p z oc r 3 - 7 ' ±0 - 40 with a correlation coefficient 
of 0.63. 

Next, we explore the possible correlations among E- lso , 
E PjZ , and To, or among Ei so , E PjZ , and Fo, using the 
stepwise regression analysis method. Our model is 
y T (x 1 ,^ 2 ) = a + bxi +cx 2 , where y, x\, and X2 stand for 
the observables (logarithmic). We measure the tightness 
of y r (x 1 , 2 : 2 ) with the dispersion and linear coefficient of 
the pair correlation between y T (x 1 , £ 2 ) and y. Our results 
are also reported in Table 2. Interestingly, much tighter 
correlations are found in our multiple variable regression 
analysis from the variable set {Ej so , E PiZ , To} than the 
pair correlations involving L- lso . The dispersions of the 
Ej so — E p z and L lso — Fo pair correlations are 0.56 and 
0.67, and their linear coefficients are 0.89 and 0.84, re¬ 
spectively. The dispersion of the relation L r p iso (E p ^ z ,T 0 ) 
is reduced to 0.33 and the linear coefficient increases to 
0.96. Figure 2 shows y r (x 1 , £ 2 ) as a function of y for the 
multi-variable correlations derived from the variable set 
{Ajso, E pz , To}. 

The multi-variable correlations derived from the vari¬ 
able set {Ejso, E PiZ , To} are less significantly improved 
over the pair correlations involving Ei so - The dispersions 
of the Ejso ~ E PiZ and E- lso — To pair correlations are 
0.57 and 0.63, and their linear coefficients are 0.83 and 
0.77, respectively. The dispersion of the E; r so (E P;Z , r 0 ) 
relation is reduced to 0.41 and the linear coefficient is 

6 Ghirlanda et al. (2012) used a different method to estimate the 
To values. Their method applies the Blandford-Mckee (1976) self¬ 
similar deceleration solution and extrapolates backwards to derive 

r 0 . 


3 


slightly increased to 0.88. We also show Ef so (E PtZ ,ro) as 
a function of Ei so in Figure 2. However, the dependence 
between E PjZ and To in the relations of E pz (Ei ao , To) 
and r 0 (-Eis O , E PtZ ) is negligible. This could be due to 
the large dispersion of the E p ^ z — To dependence in the 
current sample. 

Our analysis shows that the dispersions of the L; so — 
E P}Z and E- lso — E p , z pair correlations are comparable. 
Nava et al. (2012) reported that the L lso — E p>z cor¬ 
relation has slightly larger scatter than the Ei so — E PyZ 
correlation. This would be due to the variation of small 
sub-samples used for analysis. Note that E- lso is a time- 
integrated quantity. E p usually evolves significantly 
within a GRB and even within a GRB pulse (e.g., Lu 
et al. 2012). The E p of time-integrated spectrum is 
usually consistent with that derived from the spectrum 
around the peak time (Lyu et al. 2015). Therefore, the 
time-integrated effect of -E) so would result in significant 
scatter to the Ei SO — E PtZ correlation. It was suggested 
that the Li so — E PtZ correlation may be the basis of the 
time-integrated E- lso — E PtZ correlation (Lu et al. 2012; 
Lyu et al. 20150. The 3-parameter correlations involv¬ 
ing To improve significantly when using Li so . This may 
be due to the fact that L- lso has a double dependence on 
To (both photon energy and time), but £) so only has a 
dependence on To for photon energy. 

We focus on the {Li so , E PjZ , T 0 } correlation in our fol¬ 
lowing analysis. This correlation takes three forms, but 
two forms are of theoretical interest in the GRB physics 
(e.g., Zhang & Meszaros 2002) and GRB cosmology (e.g., 
Dai et al. 2004). 

First, in terms of L; so , it can be expressed as 
iiso,52 = io- 6 - 38±0 - 35 (F;p, z /keV) 1 - 34±0 - 14 rJ- 32±019 , (1) 


with a correlation coefficient of 0.96 and a dispersion of 
S = 0.33 dex. By adding r 0 , this relation is tighter than 
the Li so — E P}Z (Yonetoku) relation and -E) so — E P:Z (Am- 
ati) relation. It suggests that the relatively large disper¬ 
sion in the Yonetoku-relation and Amati-relation is likely 
due to the lack of consideration of the role of r 0 , a key 
parameter for the relativistic outflows of GRBs. 

The second format of the correlation is expressed in 
terms of E p ^ z , which reads 


77 . -I n 3.71±0.38i <7 „ 7-0.55it0.06-p—0.50±0.17 / 0 \ 

E PiZ = 10 keV X A; 52 1 0 , (2) 


with a correlation coefficient of 0.92 and a dispersion of 
S = 0.20. This format more directly carries a physical 
meaning, which we discuss in the next section. 

Our analysis is based on the GRBs whose To are mea¬ 
sured with the deceleration peak time in the afterglow 
lightcurves, except for GRB 060218. This method is the 
most robust one to estimate To- For a self-consistency 
check, we also introduce other GRBs whose To are de¬ 
rived from other methods. For example, Tang et al. 
(2015) systematically searched for the high energy spec¬ 
tral cutoffs from the joint LAT/GBM observations of 
GRBs and estimated r 0 for 9 GRBs using the opac¬ 
ity argument. They are consistent with the derived 


7 Since GRB lightcurves tend to peak at different time in dif¬ 
ferent energy bands, and, moreover, the time on which L; so is cal¬ 
culated is often not homogeneous in the cosmological rest-frame. 
This would make extra scatters to the Li so — E P:z correlation. 


Tiso(A PjZ , To) correlation, as shown in Figure 2. 

With the measured To for the GRBs in our sample, 
we further investigate the possible correlations among 
L iso (or E iso ), E p . and To in the jet-comoving frame 
(e.g., Ghirlanda et al. 2013), where L iso = Li so /Fp, 
E' iso = Ei so /Tq, e' p = -E p , z /r 0 . Our regression anal¬ 
ysis results are reported in Table 3. It is found that 
the pair correlations are much weaker than those among 
Liso (or Eiso), E p>z , and To and their chance probabil¬ 
ities are larger than 10 -4 . These results suggest that 
the observed pair correlations are likely governed by the 
Doppler boosting effect. The three-parameter regressions 
show that the L iso — E p ~T o relation still exists, as shown 
in Figure 3. The difference between the L iso — E p — T 0 
and the Li so — E p ^ z — To relations is the index of To. This 
is reasonable since L iso and E p are corrected by the To 
factor. 

Note that GRB 060218 is peculiar with extremely 
low L^o, E p , and To in comparison with the typical 
GRBs in our sampl<0. The inclusion of this event sig¬ 
nificantly expands the dynamical range of the correla¬ 
tions discussed above. Excluding this event, we find 
that the Lj so — E PyZ — r 0 relation still holds, showing as 

L iso , 52 = 10- 5 - 36±0 ' 49 (L;p iZ /keV) 1 - 23±0 - 14 rJ- 00±0 ' 21 and 
£p iZ = 10 3 57±o - 39 keV x T?o 59 5 t 0 ' 07r o°' 45±0 ' 17 - Within 
error bars, the indices are consistent with that derived 
from our GRB 060218-included sample. This indicates 
that the Lj so —E PjZ —To relation is not from broad ranges 
being due to this peculiar event. We also make correla¬ 
tion analysis among L iso (or E iso ), E p , and Fo by ex¬ 
cluding GRB 060218 from our sample. The results are 
reported in Table 4. One can find that the L iso — E p — r 0 
relation still exists (see Figure 3), but the indices of r 0 
are getting smaller and have a larger uncertainty than 
that in the Lj so — E p ^ z — To relation. It is also interesting 
that an L iso — E p relation is found for typical GRBs, i.e., 
logiiso .46 = (1-19 ± 0.08) + (1.18 ± 0.13) x log(Lp/keV) 
with a Spearman linear correlation coefficient of 0.85 and 
chance probability p < 10 -4 , but GRB 060218 deviates 
from this relation at a 3 a confidence level (Figure 3). 
This may hint that the radiation physics of this event 
may be different from typical GRBs (e.g., Campana et 
al. 2006; Wang et al. 2007). 

4. PHYSICAL IMPLICATIONS 

Different GRB prompt emission models have different 
predictions of the rest-frame peak energy A PjZ . Zhang & 
Meszaros (2002) summarized these predictions in their 
Table 1. In general, E PiZ is not only a function of the 
outflow luminosity L, which may be proportional to the 
observed luminosity L Pi i so , but also depends on Fo- The 
£ PiZ predictions of different models have different indices 
on L and To- With the available correlation (Eq. [2j), one 
can place strong constraints on various models. 

8 It was suggested that this event may be of a distinct population 
different from the typical GRBs (e.g., Soderberg et al. 2006, Liang 
et al. 2007), but it follows the E PtZ —Ei so relation and the spectral 
lag-L; HO relation of typical GRBs (Amati et al. 2007; Liang et al. 
2006). 
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An immediate inference from Table 1 of Zhang & 
Meszaros (2002) is that the predicted E p , z of the simplest 
baryonic photosphere models (called “innermost model” 
in the paper) and the external shock models all have 
wrong dependence on both L and T 0 . The external shock 
model is now essentially ruled out with the Swift early 
X-ray afterglow data that show a steep decay phase (see 
Zhang et al. 2006 for a detailed discussion). For the pho¬ 
tosphere model, one has E PjZ oc Tq (central engine tem¬ 
perature) oc L 1 / 4 if the photosphere radius R p h is below 
the coasting radius R c . There is no To-dependence in 
this regime, which is not supported by the data. When 
Rph ^ Rc: one has E pz oc T p ^ — Tq^R^/R f) ^ oc 

LT{L W /Tt)-^ oc L~^T^ (Zhang & Meszaros 
2002), where L w is the wind luminosity. This apparently 
conflicts with the observations. Noticing that the photo¬ 
sphere luminosity scales as L ~ L p h ~ L w (R p h/R c )~ 2 ^ 3 , 
one gets L w oc L 3 F(/ 8 . This gives 

E p , z oc T p h oc L~ 5 / 4 T 6 0 , (3) 

which violates the observations badlj0. For dissipative 
photospheres, E pz may be different from the photo¬ 
sphere temperature estimated above. However, there is 
no straightforward physics to drive E p z to satisfy Eq. ©■ 
We therefore conclude that the data do not favor the pho¬ 
tosphere models of GRB prompt emission. 

A favorable model is synchrotron radiation in an in¬ 
ternal emission region (internal shocks, Meszaros et al. 
1994; or an internal magnetic dissipation site, Zhang & 
Yan 2011). Within this model, one can generally write 
(Zhang & Meszaros 2002) 

£ p , z oc L 1/2 iT\ (4) 

where R is the radius of the emission region from the 
central engine. One can see that the 1/2 index for L 
match Eq. ©, which suggests that synchrotron radiation 
is likely the right emission mechanism. If the outflow 
is not magnetized, the emission radius is the internal 
shock radius, i.e. R ~ T^cSt oc Tq. From Eq.JU) one 
gets E PjZ oc L 1 / 2 r- 2 . The index for r 0 is too steep (-2) 
compared with the data (—0.50±0.17) (Eq©. This sug¬ 
gests that there must be another factor that moderates 
R to have a shallower dependence on To. If the outflow 
is Poynting-flux-dominated, so that synchrotron emis¬ 
sion comes from the ICMART (internal-collision-induced 
magnetic reconnection and turbulence) region, then the 
predicted E PjZ has an extra dependence on er, i.e. (Eq. 
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9 Fan et al. (2012) argued that the photosphere model can in¬ 
terpret four two-parameter correlations. However, the model can¬ 
not interpret the three-parameter correlation discovered in our pa¬ 
per. Their Eq.(3) is essentially E p , z oc L _ 1 / 4 r 2 Y^/ 2 . Noticing 
Yb = L/L w , one again derives Eq. 


58 of Zhang & Yan 2011) 

E PjZ oc L 1/2 R _1 cr 2 . (5) 

This extra dependence can make it possible to interpret 
the observed correlation. In particular, the data demand 
R~ 4 a 2 oc r-°- 50±0 - 17 _ Zhang & Yan (2011) argued that 
a larger a tends to increase R , since more collisions are 
required to finally reach the critical point for ICMART 
discharge, so that R _1 a 2 should have a shallower To de¬ 
pendence than r</ 2 . Even though more detailed numeri¬ 
cal modeling is needed to reveal the nature of our corre¬ 
lation ©, this qualitative picture seems to be consistent 
with the data. 

5. CONCLUSIONS AND DISCUSSION 

We presented a multiple linear regression analysis to 
key observables of the GRB outflows, i.e., L- lso (E - lsD ), 
Ep^, and To- The analysis of pair correlations among 
these observables well confirms several previously re¬ 
ported correlations. Most interestingly, we find a new 
tight correlation among Li so , E P}Z , and To from our 
multiple regression analysis. We argue that this tight 
iiso(E PjZ , To) correlation is more physical than the Li so — 
E P}Z and Li so — Fo correlations, and it may be directly 
related to both radiation physics and bulk motion of the 
outflows. 

The tight E PyZ (Li SO: r 0 ) correlation sheds light into the 
origin of GRB prompt emission. We show that the pho¬ 
tosphere models have difficulties to account for the cor¬ 
relation, and synchrotron radiation is likely the radia¬ 
tion mechanism of GRB prompt emission. This is con¬ 
sistent with other independent arguments in favor of syn¬ 
chrotron radiation (Zhang et al. 2013; Wang et al. 2014; 
Uhm & Zhang 2014; Zhang et al. 2015). Within the syn¬ 
chrotron model, the internal shock model predicts a too 
steep dependence on To, and hence, not favored. The IC¬ 
MART magnetic dissipation model (Zhang & Yan 2011) 
has the general merit to account for the correlation, even 
though no quantitative proof to the correlation is avail¬ 
able. 
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TABLE 1 

The REDSHIFT (z), ISOTROPIC LUMINOSITY (Li so ) AND ENERGY Ei so AS well as the peak 
ENERGY OF THE vfv SPECTRUM ( E p ) OF PROMPT GAMMA-RAYS, AND THE INITIAL LORENTZ 
FACTOR OF GRB FIREBALL (To) FOR THE GRBS IN OUR SAMPLE. 


GRB 

z 

L iso (10 52 erg/s) 

Eiso(10 52 erg) 

K (keV) 

r 0 

990123 

1.6 

27.5 ± 1.2 

356 ± 7 

1334+H 

600 ± 80 

050820A 

2.615 

o OO+0.30 

O.zz_o 45 

15-91+ 

889+H 

282^1® 

050922C 

2.198 

6-6++ 

5-44±g;2S 

629+1 

274 

060210 

3.91 

7-40+9® 

trc 9+I8.4 

709+1964 
' oz -172 

264 ±4 

060218 

0.0331 

4.27+1+ x 10~ 6 

93.2+! x 10 -5 

5.1 ± 0.3 

2.3 

060418 

1.489 

1 9 +on 
± ’ y — 0.15 

14.3 ±0.4 

572.5 

263+ 3 

060605 

3.78 

0.95 ±0.15 

2.83 ±0.45 

490 ±251 

197+° 

060607A 

3.082 

2.0 ±0.27 

9+ 

575 ± 200 

296+7 

060904B 

0.703 

0.074 ±0.014 

0.364 ± 0.074 

135 ±41 

108 ± 10 

061007 

1.261 

14.3+J 

95.4+1 

902+1 

436 ±3 

061121 

1.314 

14.1 ±0.15 

26.1 ±3 

1289 ± 153 

175 ± 2 

070110 

2.352 

0.451 ± 0.075 

5.5 ± 1.5 

370 ±170 

127 ±4 

070208 

1.165 

0.093 

0.295+1! 

66+7 

H5+0 

070419A 

0.97 

0.0098 

1 74+ 1 - 70 

97+16 

-19 

91+ 1 

071010B 

0.947 

0 'S'S'S”*” 0 
u.ooo_ 0 0 26 

9 ^?+0.19 

Z,D -0.14 

101 ± 13 

209 ±4 

071112C 

0.822 

1.047 

61.66 ±36.91 

422+7 

244 

080129 

4.394 

2.69 

7 

1349+ 2643 

j.o^y_ 80 9 

65 

080319C 

1.95 

9.5 ±0.12 

15 ± 0.79 

1752 ± 505 

228 ± 5 

080710 

0.845 

0.079 

0 - 8 ±g;| 

ouu -200 

63+ 

080810 

3.35 

9.56 ±0.83 

40.5 ± 2.9 

1364 ± 320 

409 ± 34 

081008 

1.967 

0.55 ±0.01 

a a -, +1.62 

6 - 61 ll .20 

267 ±ll 5 

250 

081109A 

0.98 

0.20 ±0.03 

4 07+ 1 - 91 
*- u ' -1.70 

208+!!® 

68 

081203A 

2.1 

2.81 ±0.19 

35 ±3 

1541 ± 757 

219+7 

090102 

1.547 

K oq+0.84 

o.oo_082 

21.5 ± 2 

1149+11! 

61 

090424 

0.544 

2-07++ 

4 q^+0.16 

*- oo -0.15 

273 ± 5 

300 ± 79 

090812 

2.452 

ioi?:§ 

42.1 ±5.5 

ro _ 549 

501 

091024 

1.092 

1.00 ± 0.22 

28 ±3 

794 ± 231 

69 

091029 

2.752 

1.72 ±0.10 

q c+0.86 
y -°-0.48 

230 ± 66 

221 

100621A 

0.542 

0.316 ±0.024 

4.37 ± 0.5 

146 ± 23.1 

52 

100728B 

2.106 

1.86 ± 0.12 

3 ±0.3 

404 ± 29 

373 

100906A 

1.727 

2.45 ± 0.09 

33.4 ± 3 

158 ± 16 

369 

110205A 

2.22 

2.50 ±0.34 

56 ±6 

715 ± 239 

177 

110213A 

1.46 

2.09 ±0.06 

6.4 ±0.6 

241±13 

223 

121217A 

3.1 

3.51 ±0.54 

62 

754 ± 230 

247 
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TABLE 2 

Results of Our Linear Regression Analysis for L lso , £ PjZ , and To in the observer frame with all GRBs in our sample, in 
which r is the Spearman correlation coefficient, p is the chance probability, and 5 is the dispersion of the reported 

relations. 


Relations 

Expressions 

r 

P 

5 

^-'iso (-^iso ) 

log Ai SOi 52 = (—0.99 ± 0.13) + (1.12 ± 0.10) x log £i ao ,52 

0.90 

< 10" 4 

0.54 

-^isoC-^Pjz) 

logL iso ,52 = -(5.10 ± 0.47) + (1.97 ± 0.18) x log(£ p , z /keV) 

0.89 

< 10~ 4 

0.56 

Liso (7)) 

log L iso ,52 = -(5.35 ± 0.63) + (2.42 ± 0.28) x logTo 

0.84 

< 10“ 4 

0.67 

-E'iso(-E'p,z) 

log Liso,52 = -(2.95 ± 0.48) + (1.49 ± 0.18) x log(£ p , z /keV) 

0.83 

< 10" 4 

0.57 

E iao (r 0 ) 

log Liso,52 = -(3.06 ±0.59) + (1.78 ±0.26) x logTo 

0.77 

< 10" 4 

0.63 

£ P ,z(r 0 ) 

log(£ p , z /keV) = (0.77 ±0.40) ± (0.82 ±0.18) x logT 0 

0.63 

< 10" 4 

0.43 

q a 0 (Ep„ r o) 

iogLf go 52 = -(6.38 ± 0.35) ± (1.34 ± 0.14) x log(£ p , z /keV) ± (1.32 ± 0.19) x logT 0 

0.96 

< 10“ 4 

0.33 

^z(Liso,r 0 ) 
rs(Li S o, Ep,z ) 

log(E p z /keV) = (3.71 ± 0.38) ± (0.55 ± 0.06) x logL ia0iS2 - (0.50 ± 0.17) x logT 0 

0.92 

< 10- 4 

0.20 

log Tg = (3.33 ± 0.38) ± (0.46 ± 0.07) x logL ia0i5 2 - (0.43 ± 0.15) x log(£ p , z /keV) 

0.88 

< 10“ 4 

0.18 

Ef so (Ep,„r 0 ) 

iogL/f a o,52 = -(3.81 ± 0.47) ± (0.98 ± 0.19) x log(£ p , z /keV) ± (0.97 ± 0.25) x logT 0 

0.88 

< 10~ 4 

0.41 

E p,z( E iso,r 0 ) 

log(££ z /keV) = (2.17 ±0.41) ± (0.46 ± 0.09) x log£ ia0i52 ± (0.00 ± 0.21) x logT 0 

0.82 

< 10~ 4 

0.26 

r^(-Eiso, -^p,z) 

log Tq = (1.92 ± 0.34) ± (0.33 ± 0.09) x log£ ia0i5 2 ± (0.00 ± 0.15) x log(£ p , z /keV) 

0.77 

< 10“ 4 

0.21 


,TABLE 3 

Results of Our Linear Regression Analysis for L iao , E p , and To in the jet co-moving frame with all GRBs in our sample, 

IN WHICH r IS THE SPEARMAN CORRELATION COEFFICIENT, p IS THE CHANCE PROBABILITY, AND <5 IS THE DISPERSION OF THE REPORTED 

RELATIONS. 


Relations 

Expressions 

r 

P 

<5* 

L iso( E iso) 

logL; s0i46 = (1.08 ± 0.11) ± (0.73 ± 0.11) x logE^ 

0.75 

< 10" 4 

0.46 

lL(K) 

lug Lj ao 4 g = -(0.31 ± 0.27) ± (0.89 ± 0.47) x log(£^/keV) 

0.32 

0.07 

- 

C(ro) 

l°gL( S o,4 6 = (0.65 ± 0.63) ± (0.42 ± 0.28) x logTo 

0.26 

0.14 

- 

Ko(K) 

log -E; so , 5 o = (0.38 ± 0.14) ± (0.82 ± 0.25) x log(E p /keV) 

0.50 

2.36 x 10“ 3 

- 

<o( r 0) 

log<0,50 = -( L06 *0-60) ± (0.78 ± 0.26) x logT 0 

0.47 

5.39 x 10“ 3 

- 

E'Fo) 

log(L p /’keV) = (0.77 ± 0.40) - (0.18 ± 0.18) x logT 0 

-0.17 

0.33 

- 

^iso (-Ep.ro) 

l°gLjg 0 jg = -(0.38 ± 0.35) ± (1.34 ± 0.14) x log(L p /keV) ± (0.66 ± 0.15) x logTo 

0.87 

< 10" 4 

0.30 

< r (E; ao ,r 0 ) 

log(Lp7keV) = (0.42 ± 0.22) ± (0.55 ± 0.06) x logL; so 46 - (0.41 ± 0.10) x logT 0 

0.86 

< 10- 4 

0.19 

r t 0 (l[ so ,e' p ) 

log Tg = (1.63 ± 0.17) ± (0.59 ±0.13) x logL| so 46 - (0.89 ± 0.21) x log(Lp/keV) 

0.64 

< 10" 4 

0.21 

-EisoC-Ep.To) 

logLi ao ,5o = -(1-81 ± 0.47) ± (0.98 ± 0.19) x log(Ep/keV) ± (0.96 ± 0.20) x logTo 

0.76 

< 10" 4 

0.35 

E P (-ELo. r o) 

log(l?L r /keV) = (1.25 ± 0.32) ± (0.46 ± 0.09) x log£^50 ± (0.54 ± 0.15) x logT 0 

0.68 

< 10- 4 

0.22 

r&(-E l so ,E' p ) 

logTg = (2.12 ± 0.08) ± (0.44 ± 0.09) x log £^0, 50 “ (°- 53 ± °- 15 ) x logl-E^/keV) 

0.66 

< 10" 4 

0.21 


Reported only for the relation with p < 10 4 


( , TABLE 4 

Results of our linear regression analysis for L iso , £ p and To in the jet co-moving frame by excluding GRB 060218 from 

OUR SAMPLE , IN WHICH r IS THE SPEARMAN CORRELATION COEFFICIENT, p IS THE CHANCE PROBABILITY, AND 5 IS THE DISPERSION OF 

THE REPORTED RELATIONS. 


Relations 

Expressions 

r 

P 

5* 

Eiso(-Eiso) 

logL; so 46 = (1.13 ± 0.13) ± (0.67 ± 0.14) x log£ ia0i5 o 

0.66 

< 10" 4 

0.46 

lUK) 

lo g Ei ao ,46 = (no ± 0.08) ± (1.18 ± 0.13) x log(£p/keV) 

0.85 

< 10" 4 

0.41 

L iso( r 0) 

logL| so 4 g = (2.45 ± 0.85) - (0.36 ± 0.37) x logTo 

-0.17 

0.34 

- 

Ko(K) 

log-EL,so = (°- 46 ± °- 12 ) + (°- 77 ± 0-20) x log(£p/keV) 

0.57 

6.00 x 10" 4 

- 

KoPo) 

log £i' so ,5o = (°- 22 ± 0-84) ± (0.23 ± 0.37) x logT 0 

0.11 

0.53 

- 

E'(r 0 ) 

log(£ p /keV) = (1.47 ± 0.59) - (0.48 ± 0.26) x logT 0 

-0.32 

0.07 

- 

W ao (-Ep.ro) 

logL^ jg = (0.64 ± 0.49) ± (1.23 ± 0.14) x log(£p/keV) ± (0.23 ± 0.21) x logTo 

0.86 

< 10" 4 

0.27 

< r (E; ao ,r 0 ) 

log(£p7keV) = (0.02 ± 0.35) ± (0.59 ± 0.07) x logL( so 46 - (0.27 ± 0.14) x logT 0 

0.87 

< 10" 4 

0.19 


log Tg = (2.15 ±0.19) ± (0.17 ±0.15) x logL^g - (0.42 ± 0.21) x log(£p/keV) 

0.37 

0.03 

- 

< r o(-E;,r 0 ) 

l°g£i' s r o , 50 = -(1.13 ± 0.72) ± (0.91 ± 0.20) x log(£p/keV) ± (0.67 ± 0.30) x logTo 

0.66 

< 10" 4 

0.30 

K 1 (-ELo.ro) 

log(£7/keV) = (1.37 ± 0.46) ± (0.45 ± 0.10) x log£i gOi 50 - (0.58 ± 0.20) x logT 0 

0.68 

< 10" 4 

0.22 

r liKo’K) 

logTg = (2.26 ± 0.07) ± (0.21 ± 0.09) x log £(^50 - (0.37 ± 0.13) x log^/keV) 

0.48 

4.91 x 10“ 3 

— 


Reported only for the relation with p < 10 4 
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Fig. 2. — Three-parameter correlations derived from our multiple regression analysis. The parameters based on the correlations as marked 
in each panel are compared against the observed parameters. The best fit lines together with their 2 a dispersion regions are shown with 
solid and dashed lines, respectively. For a self-consistency check, 9 GRBs whose To are derived with the opacity argument by using the 
high energy spectral cutoffs from the joint LAT/GBM observations are also shown (blue dots, taken from Tang et al. 2015). 





log E '(keV) 


, — E p — To relations in the jet co-moving frame. The symbol style is the same as Figure 2. An L iso — E p 
relation is found for typical GRBs, i.e., logL iso 46 = (1.19d=0.08) + (1.18±0.13) x log(l£ p /keV), but GRB 060218 deviates from this relation 
at a 3 (7 confidence level. GRB 060218 shares the same L iso — E p — Tq relation with typical GRBs. 


















